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LOCALLY D-OPTIMAL DESIGNS BASED ON A CLASS OF 
COMPOSED MODELS RESULTED FROM BLENDING 
! EMAX AND ONE-COMPARTMENT MODELS 1 

O 

CN ■ By X. Fang and A. S. Hedayat 

■ University of Illinois, Chicago 

A class of nonlinear models combining a pharmacokinetic com- 
partmental model and a pharmacodynamic Emax model is intro- 
duced. The locally D-optimal (LD) design for a four-parameter com- 
posed model is found to be a saturated four-point uniform LD design 
with the two boundary points of the design space in the LD design 
' support. For a five-parameter composed model, a sufficient condition 

for the LD design to require the minimum number of sampling time 
points is derived. Robust LD designs are also investigated for both 
models. It is found that an LD design with k parameters is equivalent 
to an LD design with k — 1 parameters if the linear parameter in the 
two composed models is a nuisance parameter. Assorted examples of 
LD designs are presented. 



qq ■ 1. Introduction. A class of models is constructed by plugging a pharma- 

\ cokinetic (PK) compartmental model into a pharmacodynamic (PD) Emax 

model. Under this class of models, only one measurement is required per 
study subject rather than multiple measurements and both PK and PD pa- 



(N 



rameters can be estimated by a single experimental setup. Other advantages 
■ of this approach will be listed shortly after the related PK and PD models 

O ! are introduced. 

A basic PD Emax model can be expressed as E([D]) = Eq + E max [D]/ 
(ED$q + [D]), where E(-) is the drug effect such as reduction in low-density 
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cholesterol and [D] is the concentration of free drug in the environs of the 
drug receptor. The Emax model contains three PD parameters. ED^q is the 
drug concentration showing 50% of the maximum drug effect, E mSuX is the 
maximum drug effect and Eq is the baseline effect. For the concept of Emax 
model, the reader is referred to Ritschel [22]. Various applications of Emax 
models have been discussed by, but are not limited to, Graves et al. [7], 
Demana et al. [4], Angus et al. [1], Staab et al. [26] and Hedayat, Yan and 
Pezzuto [11, 12]. 

Since a PK compartmental model describes drug concentration across 
time, the [D] can be approximated by an appropriate PK model. After 
replacing [D] in the Emax model by an open one-compartment model with 
IV bolus input and first-order elimination results in the composed Emax- 
PK1 model 

(1-1) Y >i = A + p^M^Td + e <!< 

where Ytj is the effect of the drug at time t observed on subject j, (3q is the 
minimum residual effect, (3\ is the maximum drug effect, /?2 is the ED§$ and 
(3?, is the total elimination rate. It is assumed that the errors across time and 
subjects are i.i. d. N(0,a 2 ). Throughout this paper, D is the administered 
dose in the unit of drug amount per unit body weight. For the concept 
and essential results of the one-compartment model related to IV bolus, 
the reader is referred to Rowland and Tozer [24], Landaw [16], Dette and 
Neugebauer [5], Han and Chaloner [8] and Hedayat, Zhong and Nie [10]. 

If [D] is replaced by an open one-compartment model with the first-order 
input and the first-order elimination, the resulting composed Emax-PK2 
model becomes 

01 ■ D(e~ l32t - e~P 3t ) 

(L2) Ytj =P0 + Ml-fr^ + Die-M-e-M) + % ' 

where (5q is the baseline effect, (3\ is £ max , /?2 is the absorption rate, /?3 is the 
total elimination rate and (5^ is the ED§q. The assumption about the error 
terms is the same as that in model (1.1). For the concept and various results 
related to this one-compartment model, the reader is referred to Rodda, 
Sampson and Smith [23], Davidian and Gallant [3], Mallet [18], Mandema, 
Verotta and Sheiner [19], Verme et al. [27], Lindsey et al. [17], Landaw [16], 
Atkinson et al. [2] and Wakefield [28]. 

Other advantages of the Emax-PK models include (1) The PK parameters 
can be estimated without drawing blood samples. (2) Drug effect becomes a 
function of time rather than a function of drug concentration which itself is a 
function of time. As a result, the drug effect is predictable across time, such 
as blood pressure is being reduced across time. (3) Sampling at various time 
points is controllable before taking samples, whereas the drug concentration 
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in an Emax model alone is unknown before sampling. Therefore, one can 
obtain better estimates of PK and PD parameters by implementing suitable 
designs. (4) The residual effect of a drug can be estimated at any time point 
through the Emax-PKl model. (5) The baseline effect of a drug can be 
estimated through the Emax-PK2 model. 

In this paper, locally D-optimal (LD) designs and robust LD designs for 
the preceding two Emax-PK models are studied for the purpose of estimat- 
ing all population PK parameters. These PK parameters are considered fixed 
effects although they likely differ between subjects. In Section 2, the equiv- 
alence theory for nonlinear models is briefly reviewed. Section 3 contains 
the main results about the number of sampling time points in an LD design 
support for the Emax-PK models. Examples of LD designs for both models 
are provided in Section 4. Some assorted robust LD designs are investigated 
in Section 5. Discussion and conclusion are summarized in Section 6. 

2. Preliminaries. The equivalence theorem for designs based on linear 
models was first introduced and developed by Kiefer and Wolfowiz [15]. 
White [29] extended the equivalence theorem to locally optimal designs 
based on nonlinear models. In a nonlinear model, the Fisher information 
matrix depends on unknown parameters. As a result, there is no global op- 
timization for all values of model parameters and the optimality of a design 
can be evaluated only for special values of model parameters, called nomi- 
nal or postulated values. For details, the reader is referred to Silvey [25] and 
Hedayat [13]. 

Without loss of generality, let us consider the nonlinear model Ytj = 
rj(t,P) + Etj with Em's being i.i.d. N(0,a 2 ). The normalized Fisher infor- 
mation matrix for the entire model parameters based on the design measure 
£ = {(t\,pi), . . . , (tx,pK)} can be expressed as 

A design measure £ = {(t±, pi ),■■■, (tK,PK)} is a description of sampling 
time points (t\, . . . ,tx) where Ytj will be measured for the jth subject at 
time ti. Associated with U is the mass pi such that < pi < 1, and J2n = 1- 
In design £, pi, ■ ■ ■ ,pk represent the proportion of the number of subjects 
studied taken at time ti,...,tx, respectively. For a total sample of size n, 
rii = npi is the number of subjects to be studied at ti. For the purpose of 
obtaining the most precise estimators of the entire model parameters, one 
needs to identify a design measure whose related Fisher information matrix 
is nonsingular and whose determinant is maximized in the class of competing 
designs. This is because M~ 1 (^, f3) is proportional to the asymptotic variance 
and covariance matrix of the MLE of the model parameters. An LD design 
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for a given set of parameters has the maximum determinant for its Fisher 
information matrix based on the postulated values for these parameters. 

When the number of PK parameters of interest is s < k in a /c-parameter 
nonlinear model, the asymptotic generalized variance of these s estimators 
can be expressed as ?i _1 [M(£,/?)/M 22 (£,/3)] _1 , where 

M 21 (£,/3) M 22 {iJ)j 

is the partitioned form of M(£,/3) with Mil (£,/?) being the associated s x 
s information matrix of the s parameters under the design £. Following 
White [29], the design £* in a class of competing designs, T>, is said to be an 
LD S design if det M(£* , 0)/ det M 22 (£* , 0) = max^p [det M(£, det M 22 
White [29] established that the design £* is an LD S design if and only if 
su PteT d(t, CJ)=s, where d(t, £, (3) = tr{/(t, /?)M _1 (£, /?)} - tr{/ 22 (i, /3) x 
M 22 1 (C,/3)} when s < /c and d(t,£,/3) = tr{/(i,/3)M _1 (£,/?)} when s = fc and 
T = [0, oo) is the set of all sampling time points with the understanding that 
the extreme right time point is a very large practitioner-selected time point. 
Here, I(t,/3) is the information matrix for the design with the support point 
of t only and it is partitioned as that of M(£,/3). The quantity d(t,£,f3) can 
be interpreted as the variance of the estimated response at time t when s = k 
and is the variance of the estimated response after eliminating the effects of 
the k — s nuisance parameters when s <k. 

Before searching for optimal t{ and pi, the required number of design 
points in an LD design support is investigated first for both theoretical and 
practical interests. Knowing the required number of design points in advance 
would significantly reduce the task of searching for an LD design. 

3. Main results. In this section, the support size of an LD design for the 
Emax-PKl model under the normality assumption is investigated. Since the 

corresponding induced design space, C = {a" 1 { ^%$ ) '-t G T}i is a bounded 

subset in M fc and the determinant of the Fisher information matrix is a 
continuous function on C, thus an LD design for this setup must exist. In 
what follows, by a saturated design it is meant a design in which the number 
of design time points is equal to the number of model parameters. 

Theorem 1. Under model (1.1): E(Y t ) = j3 + n%F +D with Var(y t ) = 

a 2 , where t £ [0,u], an LD design is a saturated D-optimal design with 
time and u in its support. 

Proof. Let M (£*,/?) be the Fisher information matrix for f3 = (J3q, /3±, /3 2 , 
(3 3 ) T based on an LD design £* and define fo(t) =ti(I(t,(3)M~ 1 (^*,(3)) + c, 
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D{f3 2 e^ + DY l 
V -faDfateM ((3 2 e^ + D)~ 2 j 



D{f3 2 e^ t + DY l 
f3 l De^ t {f3 2 e^ t + D)~ 2 
\-f3 1 Df3 2 te^ t ((3 2 e^ t + D)~ 2 J 

It can be shown that fo(t) has the same number of zeros as f(t) = a^e^ 1 + 
e 3A * (a 31 t + a 30 ) + e 2 ^ 1 {a 22 t 2 + a 2 it + a 20 ) + e fti (aut + a 10 ) + a 00 , where ay- 
depends on and /3 fc , i G {0,1,2,3,4}, j G {0,1,2}, k € {1,2,3}. No- 
tice that a 22 > and 031 < since the cofactor Cofi4 in M is positive 
(Appendix). Also, by the Remark in the Appendix, it can be shown that 
an > 0. Now in order to prove that f(t) has at most six zeros, the proper- 
ties of various derivatives of f(t) will be explored. 
By differentiating f(t) with respect to t, one has 

f'(t) = 4/3 3 a 4 oe 4/33t + e 3/33t (3/3 3 a 3 it + 630) 

+ e 2Ai (2/3 3 a 22 t 2 + b 21 t + b 20 ) + e^ftani + ho), 

where the fr^'s are the corresponding constants. After setting f'(t) =0 and 
multiplying both sides of f'(t) = by e -4 ^ 3 *, one obtains 

f'(t) = 4/? 3 a 40 + e-^(3(3 3 a 31 t + 6 30 ) + e -^ 3t (2/3 3 a 22 t 2 + b 21 t + b 20 ) 

+ e- 3/33t (/3 3 ant + Oio). 

Notice that f'(t) and f'(t) have the same number of zeros. After differenti- 
ating f'(t), one has 



/"(*) 



-fat. 



-3Pia 31 t + c 30 ) + e- 2/33 *(-4/?|a 22 t 2 + c 21 t + c 20 ) 



+ e -3&t(_3^ Ollt + Cl0 ) 



where Cy-'s are the corresponding constants. By multiplying f"(t) by e 
one has 

/"(*) = e % *(-3/3fa 3 it + C3o) 

+ e /33i (-4/3 3 2 a 22 t 2 + c 21 t + c 20 ) + (-3%a n t + c w ). 

Differentiating f"(t) yields 



3fat 



fit) 



Jifati 



0ta 31 t + d 30 ) + e k \-A0^a 22 i l + d 21 t + d 20 ) - 3^ou, 
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where dij's are the corresponding constants. Next, it will be shown that the 
function 

g (t) = e 2 ^(-6f3la 31 t + d 30 ) + e ftt (-4/? 3 3 a 22 £ 2 + d 21 t + d 20 ) 

has at most three stationary points. By multiplying both sides of g'(t) = 
by e -2 ^ 3 *, one has 

g'(t) = -12/3la 31 t + e 30 + e-^ l {-^ 3 a 22 t 2 + e 21 t + e 20 ) = 0, 

where e^-'s are the corresponding constants. Setting g"(t) = 0, one has 

(3.1) e-^\Apla 22 t 2 + f 21 t + / 23 ) = 12/3 3 4 a 31 , 

where /ij's are the corresponding constants. Since the left-hand side of (3.1) 
has at most two stationary points and it approaches above the t axis as t 
goes to co, the most right monotone interval of the left-hand side will not 
intersect the horizontal line y = 12j3 3 a 3 \, which is below the t axis. Therefore, 
(3.1) has at most two roots. 

Now, since (3.1) has at most two roots, g'(t) will have at most three zeros. 
This implies that g(t) has at most three stationary points. Since g(t) goes to 
below the t axis as t goes to — oo and 3/3|an > 0, the most left monotone 
interval will not intersect the horizontal line y = 3/3 2 an, which is above the 

t axis. Consequently, f"'(t) has at most three zeros. This implies that f"(t) 
has at most four zeros and f'(t) has at most five zeros. As a result, f(t) has 
at most six zeros. 

Since f(t) has at most six zeros for all c's and by the equivalence theorem 
the interior points of an LD design must be locally maximum points of 
tr(/(i,/?)M -1 (£*,/3)), therefore there are at most two optimal design points 
on (0, oo). Otherwise, f(t) has more than six zeros for some c. Since the 
induced design point at time is not proportional to that at time u > 0, 
the existence of an LD design based on model (1.1) forces the two boundary 
points to be in the optimal design support. □ 

Theorem 1 demonstrates that an LD design for model (1.1) is of the 
form C = {(0, 1/4), (t 2 , 1/4), (t 3 , 1/4), (u, 1/4)} over t G [0,u\. Consequently, 
to search for an LD design for model (1.1), one only needs to find out t 2 and 
is- 

Theorem 2. An LD design for the Emax-PK2 model has minimum 
support size with time point in its support if d(t,^,P) has precisely four 
locally maximum time points when t € (0,oo). 

Proof. By Caratheodory's theorem [25], for this five-parameter nonlin- 
ear model, an LD design must have at least five time points in its support. 
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By the equivalence theorem, LD design time points must be necessarily the 
locally maximum points of d(t, £, 0) when t S (0, oo). Since d(t, £, 0) has four 
locally maximum time points and the induced design time points at t = 
and t = oo are the same, the existence of an LD design forces time point 
and the four locally maximum points to be in the design support. □ 

Theorem 2 provides a sufficient condition for an LD design for the Emax- 
PK2 model to be minimally supported when the original design space is 
[0, oo). In general, if the required number of LD sampling time points is un- 
known for a fc-parameter nonlinear model, the practitioner could search for 
the best fc-point uniform LD design first. Then the V-algorithm by Fedorov 
[6] can be applied to search for an LD design with the fc-point uniform LD 
design as the initial design. 

When 0o is a nuisance parameter, it is found that an LD design is also 
an LDfc_i design for a /c-parameter nonlinear model of the form ytj = 0o + 
fj(t,0) +£tj, where etj are i.i.d. N(0,a 2 ) and dfj(t, 0)/d0i is not free of 0i 
for i = 1, . . . ,k — 1. For a model of this form, the following result is obtained. 

Theorem 3. An LD design for the k -parameter nonlinear model is also 
an LDk_i design for parameters 0i,..., 0k— l ■ 

The proof is based on the fact that both I22 in I(t, 0) and M 2 ~ 2 1 in M(£*,0) 
are equal to 1. Theorem 3 shows that the LDfc_i design is globally optimized 
for 0o although it is only locally optimal for the other nonlinear parameters. 
Theorem 3 is applicable to both Emax-PK models discussed in this paper. 
Examples of LD^_i designs originated from LD designs for both Emax-PKl 
and Emax-PK2 models are given in Section 4. 

4. Illustrated examples. Maximizing the determinant of the Fisher in- 
formation matrix under nonlinear models requires the foreknowledge of the 
values of model parameters. In practice, a good guess can be obtained from 
a pilot experiment. For example, suppose based on the prior information 
from the pilot experiment, the minimum residual effect is equal to 0.5, the 
maximum drug effect is equal to 10, the ED$o is equal to 1 mg/kg and the 
total elimination rate is equal to 0.1 hour -1 . The search for an LD design 
generally contains two steps: (1) finding a best fc-point uniform LD design; 
(2) using the V-algorithm (Fedorov [6]) to find an LD design given the ini- 
tial design as that found in step (1). For the Emax-PK2 model, a search 
as described was performed. However, for the Emax-PKl model, the search 
was stopped at step (1) since Theorem 1 shows that an LD design for this 
model is a saturated LD design. 

For the Emax-PKl model with postulated = (0.5, 10, 1,0.1) and an ini- 
tial dose of 5 mg/kg, a four-point LD design is found to be uniform at 
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times 0, 13.48, 31.62 and 160.00 hours. Figure 1 illustrates some essential 
characteristics of d(t, £, /3) for this design. 

If the maximum sampling time is 72 or 24 hours, then the LD design can 
be found to be uniform at times 0, 13.26, 31.05 and 72 hours or at times 
0, 7.97, 17.83 and 24 hours, respectively. The d(t,^,(3) for these two designs 
are illustrated in Figures 2 and 3. 

Although the upper bound of the support size for an LD design based 
on the Emax-PK2 model is 16 by Caratheodory's theorem, one can still 
follow the two-step search described above. For the Emax-PK2 model with 
postulated (3 = (0.5,10,0.5,0.1,1), an LD design is found to be uniform at 




(a) (b) 



Fig. 1. Plots of d(t,C,0) with (a) t G [0,20] and (b) t £ [20,200], respec- 
tively. The design space is [0, oo) and the optimal sampling times are found at 
0,13.481506,31.624511,160.00199. The initial IV dose is 5 mg/kg. 




Fig. 2. Plots of d(t,C,0) with (a) t€[0,20] and (b) t 6 [20,200], respectively. 
The design space is [0, 72] hours and the optimal sampling times are found to be at 
0,13.263029,31.050686 and 72 hours. 
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times 0, 0.275, 2.999, 14.75 and 32.7 hours. The associated d(t,£,(3) for this 
design is illustrated in Figure 4. 

5. Robust designs. Applying an LD design in practice may be criticized 
for its local optimality. If the nominal values of the model parameters are not 
close to the true values, a more desirable design would be a robust design, 
which would lead to a better parameter estimation than an LD design while 
maintaining high efficiency. The relative efficiency of a design £ compared 
to an LD design £* for the postulated (5 is defined as det (£,/?)/ det (£*,/?). 

A robust index defined as |(<9(det M(£,/3))/<9/3j) -1 | for parameter i = 
0,1,... is introduced to compare the robustness of LD designs based on 
different nominal values of the model parameters. For a linear model, this 
robust index is oo, indicating that an LD design is globally optimal for all 
values of the model parameters. However, for a nonlinear model, it measures 
the inverse of the changing rate of the determinant in the neighborhood of 
the postulated (3. Therefore, it is called the locally robust index (LRI). The 
larger the value of the LRI, the more locally robust the design. 

One class of robust designs is the equally spaced uniform LD (ESULD) 
designs. An advantage of this class of designs is that it can be implemented 
very easily. It is robust for the estimation of parameters in the models that 
were described in Hedayat, Yan and Pezzuto [11, 12] and Hedayat, Zhong 
and Nie [10]. The class of ESULD designs is examined here for both Emax- 
PK1 and Emax-PK2 models. Table 1 shows that an ESULD design retains its 
robustness but loses its efficiency as the size of the design support increases. 
For practical applications, the five- or the six-point ESULD designs are 
recommended. 

Since the efficiency of the ESULD designs is very low for the Emax-PK2 
model, another class of robust designs is constructed based on an LD design. 




5 ID 15 20 ZS 



t 

Fig. 3. Plot of d(t,£* , (5) when the design space is [0,24] hours based on the Emax-PKl 
model. 



o 



Table 1 

ESULD designs for the model with the initial dose = 5 mg/kg and the nominal values /3i = 10, /?2 = 1 mg/kg, (3z = 0.1 h -1 

5 

Support size Support of £ Efficiency of £ LRI for /3i LRI for p2 LRI for j3s C 

"1 0,13.48,31.62,160.00 I 0.95901654 0.10344687 0.019187232 q 

5 0,16.176 X i,i = 1,2, 3, 4 0.476373989 2.0131589 0.19911378 0.040256987 > 

6 0,14.736 x i,i= 1,2,3,4,5 0.481038755 1.9936368 0.20439842 0.039876038 § 

7 0,13 x i,j = 1,2,..., 6 0.440686613 2.1761871 0.22903063 0.043844867 > 

8 0,11.094 x i,i = 1,2,..., 7 0.403119466 2.3789884 0.25133456 0.047580261 M 

9 0,9.647xi,i = 1,2,..., 8 0.376132417 2.549678 0.26807298 0.050989122 ffi 

10 0,8.558 x«,i = l,2,...,9 0.355042847 2.701129 0.28268401 0.054018207 H 

11 0,7.686 Xi,i = 1,2,..., 10 0.3379349 2.8378736 0.29572619 0.056762451 > 

12 0,7xi,i = l,2,...,ll 0.323818285 2.9615886 0.30723535 0.058653795 £ 

When s = 4, the design is the LD design. ^ 
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t t 

(c) (d) 

Fig. 4. Plots ofd(t,C,0) with (a) tG [0,0.5], (b) t 6 [0.5,20], (c) te [20,60] and (d) 
£ 6 [60,200], respectively. The numerical solutions of the local maximums are 0, 0.275, 
2.999, 14.75, 32.695. The initial dose is 5 mg/kg. 



The support of such a robust design includes all the design time points of 
an LD design as well as the design time points in the form of t* + r or t* — r, 
where t* is a design time point of an LD design and r is a fixed number. 
For convenience, this class of designs is referred to as equal-step expanded 
uniform LD (ESEULD) designs. For example, if an LD design is 

t* = ( *i *2 H \ 
? ~~ ^1/3 1/3 1/3 ) ' 

then a four-point ESEULD design would be 

t** = ( h h+r t 2 t 3 \ 
? v l/4 1/4 1/4 1/4 

a five-point ESEULD design would be 

t** ( h h+r t 2 h + r t 3 \ 
^ \l/5 1/5 1/5 1/5 1/5 / ' 




a six-point ESEULD design would be 

t** _ ( h h+r t 2 t 2 + r t 3 -r t 3 \ 
^ \1/G 1/6 1/6 1/6 1/6 1/6 / 

and a seven-point ESEULD design would be 

« = / ti h+r t 2 -r t 2 t 2 +r t 3 - r t 3 \ 
? ^1/7 1/7 1/7 1/7 1/7 1/7 1/7/' 

Since the relative efficiency of the design £** relative to an LD design £*, 
det M(£**, /3) I det M(£*,/3), is a function of r given £* and (5, it is impossible 
to give an explicit form of r as a function of efficiency. However, it is possible 
to plot efficiency versus step length r to find out the numerical relationship 
between r and the efficiency. Figure 5(a) and 5(b) shows such plots based 
on the Emax-PKl model and the Emax-PK2 model, respectively. The same 
postulated values of as those in Section 4 are used for illustrations and 
tables throughout this section. From this figure, it is clear that the value 
of r cannot be assigned arbitrarily for some n-point ESEULD designs. For 
example, r = 0.95 does not exist for a five-point ESEULD design based 
on the Emax-PKl model. In this figure, r = 1 hour and r = 0.2 hour are 
chosen for robustness study based on the Emax-PKl model and the Emax- 
PK2 model, respectively. The corresponding design time points are listed in 
Tables 2 and 3. The LRFs are calculated and listed in Tables 4 and 5. From 
these results, it appears that the higher the relative efficiency of a robust 
LD design, the less the robustness of the design. 



6. Conclusion and discussion. This paper introduced a class of models 
by blending a PD Emax model and a PK compartmental model and studied 
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some important features of LD, ESULD and ESEULD designs for the Emax- 
PK1 and the Emax-PK2 models in this class. 

For Emax-PKl model, an LD design is a saturated four-point uniform 
LD design with the two boundary time points of the design space in its 
support. Both time and the upper bound of the design space are the infor- 
mative time points here. This can be observed directly from the model. As 
t approaches oo, the Emax-PKl model is reduced to Y t j = (3q + etj. There- 
fore, the upper bound of the design space is an informative time point for 

Table 2 

ESEULD designs based on the Emax-PKl model with r = 1 hour 



Support size ESEULD design point 



4 




0,13.48,31.62,160.00 


5 




0,1,13.48,31.62,160.00 


6 




0,1,13.48,14.48,31.62,160.00 


7 




0, 1, 13.48, 14.48, 31.62, 32.62, 160.00 


8 




0, 1, 13.48, 14.48, 31.62, 32.62, 159, 160.00 


9 




0, 1, 12.48, 13.48, 14.48, 31.62, 32.62, 159, 160.00 


10 




0, 1, 12.48, 13.48, 14.48, 30.62, 31.62, 32.62, 159, 160.00 






Table 3 




ESEULD designs based on the Emax-PK2 model with r = 0.2 hour 


Support 


size 


ESEULD design points 


5 




0,0.275,2.999,14.75,32.695 


6 




0, 0.275, 0.475, 2.999, 14.75, 32.695 


7 




0, 0.275, 0.475, 2.999, 3.199, 14.75, 32.695 


8 




0, 0.275, 0.475, 2.999, 3.199, 14.75, 14.95, 32.695 


9 




0, 0.275, 0.475, 2.999, 3.199, 14.75, 14.95, 32.495, 32.695 


10 




0, 0.275, 0.475, 2.799, 2.999, 3.199, 14.75, 14.95, 32.495, 32.695 






Table 4 




Efficiency and robustness of ESEULD designs for the Emax-PKl model 


ESEULD design 


Efficiency of £ LRI: /3i LRI: /3 2 LRI: (3 3 


4-point 




1 0.95901654 0.10344687 0.019187232 


5-point 




0.7588 1.2638471 0.135847 0.025576467 


6-point 




0.7284 1.3165725 0.13967382 0.025674003 


7-point 




0.7846 1.2223363 0.12857595 0.021936078 


8-point 




0.9193 1.0431845 0.10970443 0.018685843 


9-point 




0.8583 1.1172873 0.1190755 0.020901414 


10-point 




0.8437 1.1366831 0.12194774 0.022945661 
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Table 5 

Efficiency and robustness of ESEULD designs for the Emax-PK2 model 



ESEULD 
design 


Efficiency of £ 


LRI: (3i 


LRI: /3 2 


LRI: f3 3 


LRI: f3 4 


5-point 


1 


0.3067854 


0.046489413 


0.0090373459 


0.043066663 


6-point 


0.7667 


0.40014335 


0.055507767 


0.011777852 


0.058844405 


7-point 


0.6995 


0.43857707 


0.059457229 


0.013001578 


0.064883822 


8-point 


0.7168 


0.42799782 


0.05768681 


0.012643369 


0.063180006 


9-point 


0.7952 


0.38580638 


0.052011575 


0.011571927 


0.057033582 


10-point 


0.6995 


0.4385802 


0.059583199 


0.013097578 


0.064855939 



parameter (3q. When t = 0, the model is reduced to Ytj = (3q + fo +D + e tj- 
Consequently, time is an informative time point for (3q as well as the ratio 
of fa /fa- 

A sufficient condition for the Emax-PK2 model to be minimally supported 
is given. Time is an informative time point here. This also can be explained 
from the Emax-PK2 model directly. As time goes to 0, the Emax-PK2 model 
is reduced to Ytj = (3q + etj ■ 

When (3q is considered as a nuisance parameter, an LD design and an 
LDfc_i design based on any of the Emax-PK models are equivalent. The 
corresponding LD design is globally optimal for the linear parameter (3$ and 
locally optimal for the other nonlinear parameters. 

Future research for the Emax-PK models could involve random effects for 
certain PK parameters since these parameters likely differ between subjects. 
The reader is referred to Mentre, Mallet and Baccar [20], Palmer and Muller 
[21], Han and Chaloner [9] for recent results in this area. 



APPENDIX: PROOFS 



Cofi4 of M(^*, ( 3)- 1 is positive and Cof 24 of M(|*,/3) _1 is negative. 

Since the Fisher information matrix of f3 = {Pq, Pi, P2, Pz) T at sampling 



time ti for an LD design £* 
Emax-PKl model is 



Vpi 



P2 



based on model (1.1) or the 



I(U,P) 



o 2 Pi 



\ 



D{p 2 e^ u + Dy 1 
PiDe^ u (p 2 e^ u + D)' 2 
\-PiDp2tie^{P 2 e^ + D)- 2 ) 

1 V 

D(p 2 e^ k + D)- 1 

-p 1 De^(p 2 e^ +D)~ 2 
PiDp 2 t i e^(p2e l33U + D)~ 2 J 
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the Cofi4 has the same sign as 
/ D 



detmi4 = — det 



Let 



i=l 
s 

E Pi i r- 



(fce^ti + D) 

-Dfce 03 ** 
(foe^+D) 2 



V^ Pl (/ 

V i= 1 



z=l 
s 

E^ 

i=l 
a 



' (/3 2 e^'i + D) 2 
' (fretoi + D) 3 
(faePaU + D)3 



i=l 

8 

E p * 

i=l 
a 

E pi 



D 2 p 2 e 203t ' 
(/3 2e /3 3 ti + D)4 

D 2 p 2 (3 2 Ue 2l33ti 



g(ti,tj,t k ) = det 



D 



D 2 



{fae^i + D) 2 



(foe®*** + Df 
tjt 



-D 2 p 1 e f3:itk \ 

D 2p2 e 2/3 3 t k 

(foePtfk + Dy 
D 2 P 2 (3 2 t k e 2 ^ 
(foefotk + D y ) 

det mi4 



mi4 



is 



\/i < j < k < s, then the determinant of 

^2i<j<k<s^2permutations of i,j t kPiPjPk9v'ii^j ■ 

Now without loss of generality, let i = 1, j = 2, k = 3; then it can be shown 
that 



E 9{h,t 2 ,t 3 ) 

permutations of 1,2,3 



D 7 f3fP 2 



(foeMi + L>) 4 ( ) 5 2 e**2 + Z?) 4 (/? 2 e^ 3 + D) 4 



At 3 



he?** 
t 2 e^ st2 



s /%*i e 2/3 3 ti 

s /3 3 *2 g2/3 3 t 2 
ol 3 3 t 3 p 2/? 3 t 3 



since 



g(h,t 2 ,t 3 ) 



D 6 pf(3 2 



(f3 2 e^ + D) 4 (p 2 e^ t2 + D) 4 (/3 2 e** 3 + D) 4 



o/5 3 <2 
o/3 3 * 3 



t 2 e^ t2 



x (/^e^ 3 * 1 +L>) 2 (/3 2 e /33 * 2 + L>)e /33 * 3 . 

Now since both {1, e^ 3 *, te^ 3 *} and {1, e^ 3 *, e 2 ^ 3 *} are Chebyshev systems (see 
Karlin and Studden [14]) and their determinants are positive, 
T. P ermutations of 1,2,3 ) *2, <3) < 0. Therefore, det 77114 = 

Epenwrtoaoni o/ i,j ,kPiPjPk9(U,tj ,t k ) < 0. This implies that the 

Cof 14 > 0. 

Next, Cof 24 < will be shown. By the same argument, it can be proved 
that 



Cof 24 = PiPjPkQ 

i<j<k<s 



ti 

h 



t k e l3ztk 



e P3U 

e Pa*i 



Cof 



14, 



16 



X. FANG AND A. S. HEDAYAT 



where 

Q = ~ {foeM* + £>) 4 (/3 2 e^ + £>) 4 (/3 2 e^ + D) A ' 

Since {l,i,e^ 3 *} is a Chebyshev system (see Karlin and Studden [14]) and 
Cofu > 0, this yields Cof 2 4 < 0. 

Remark. 

-Cofu - Cof 24 = - PiPjPkQ 

i<j<k<s 

which is positive. 
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